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Tapered Wiggler Analysis of High Grain 
Free Electron Laser Oscillators 

I. Introduction 

The development of lasers in which the active medium is a relativistic stream of free 
electrons has recently evoked much interest. The prospect of operating these continuously 
tunable, high gain devices at high efficiencies by appropriately tapering the wiggler field 
parameters has significantly increased the interest in these devices. The dramatic effect, 
of efficiency enhancement schemes in free electron lasers relies on the maintenance of the 
co-resonance between the phase velocity of the ponderomotive wave with the decreasing 
velocity of the electron beam as it gives its energy to the radiation field. Variable parameter 
wiggler free electron lasers have been previously considered for operation in the amplifier 
configuration (Refs. 1-3) and also in the micro-pulse, low gain operation in the oscillator 
configuration (Refs. 4,5). The present analysis extends the range of application to the case 
of intense, long pulsed electron beams in which the self-consistently derived space charge 
potential (arising from the axial bunching of the particles in the beat wave) is not negligible. 
Also included in this analysis are arbitrarily relativistic beams (no ultra-relativistic beam 
approximation). Thusly, the simplifying assumptions that 7 > 1 and (3. == 1 are not 
incorporated into the theory. The full relativistic effects of mildly relativistic beams, which 
are characteristic of beams from intense, long pulsed devices, are adequately considered. 
The analysis also includes the high gain operation in the oscillator configuration which is 
characterized by exponential field growth in the linear regime. 

The organization of the paper is as follows. In Sec. II, we present the model and as¬ 
sumptions used in the numerical simulation. In Sec. Ill, we derive approximate theoretical 
expressions for the nonlinear efficiency. In .Sec. IV, we present a summary of conclusions 
with specific examples from the simulation and comparison with theoretical expressions. 

II. Simulation Model 

We have conducted a fully relativistic, nonlinear simulation of the field and particle 
evolution in a variable parameter wiggler, high gain free electron laser oscillator. The 
simulation model which is kinetic in nature is coupled with Maxwell's equations to self- 

consistently include the evolution of the space charge waves. The wiggler field we consider 
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is helically polarized in the plane transverse to the beam propagation direction, and is 
ideal in that the transverse variations in the vector potential are neglected. 


A w (z) = -a u ,(;)^ exp[-i J k w {z')dz'} c+ + c.c.^, (1) 

where a u ,(r), k w (z) are the arbitrarily prescribed wiggler amplitude and wave number re¬ 
spectively and also e± - (e x ± te„)/2. and c.c. denotes complex conjugate. The radiation 
field model is also one-dimensional and is characterized by a temporally stationary ampli¬ 
tude with arbitrary spatial dependence. 

A,(z,t) = <* 4 (z)exp[*(fc 4 z - u> 4 /)] e + + c.c. (2) 

The electrostatic space charge potential is also assumed to be temporally stationary with 
arbitrary spatial dependence. 

${z,t) -- 0!(z)exp[t{ / (fc. -I- A; M .(z')|dz' u> 4 f}) + c.c. (3) 

J 0 

The vector and scalar potentials of the electromagnetic wave evolve according to Maxwell’s 
equations and the particles evolve in accordance with the kinetic distribution / (Ref. 6). 
Transverse effects are heuristically included by considering / dr f dy f 
ot Trrg, J dx J dy a, ex n with rj, equal to the electron beam radius and «i 0 = 
(2 Zft/kf) 1 / 2 is the spot size for an oscillator mode of Rayleigh range zp. Inserting the 
driving currents and charge densities into Maxwell’s equations and assuming i) a cold elec¬ 
tron beam, ii) vector and scalar potentials that vary slowly as compared to an optical 
period and iii) small radiation fields, o 4 <<r a w , one obtains the following, 


da, i ( Wb\ 2 Fi3« „ / e ' v ’ \ 

nr = 2 It) 

/ d i{k, + K, )\. i/* h y lh o 

\dz + 2 r l 2V c) (k. + k v )\ & r 


(4) 

(5) 


where we have introduced the dimensionless potentials, d.,, a w , 4>\ corresponding to 
|e|(a,,Ou,,^i)/mc 2 . The nonrelativistic beam plasma frequency is given by u>(, = 


2 






(4ir»e 2 /m ) 1 ^ 2 and F = t\/wq is the filling factor. The integration over initial positions 
has been converted to an integral over initial phases, ij> 0 , with {•••) = (1/2tt) J* w ■ ■ • dt/’o, 
and the phase is given by ^ [fc, -f k u - u > t /v z (z')\dz' + V’o- 

The evolution of the particle phase as determined by the Lorentz force equation is 
given by, 


r,« _ fgf + ^^±- k > _ 

w vj dz w 2f3 3 zl 2 dz 


k» dn", (1 4' fit)k a k w a w . 

~~2i {a ‘ r ~ "- e 


-></• 
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*.(!-#) 




+ MM + *-*( 4 : ~ *(*«. + MM 
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where prime denotes differentiation with respect to z. The oscillator system of equations 
(Eqs. (4,5)) are solved with the following methods and boundary conditions. The ordinary 
differential equations for the fields and particles are solved with a four point Adams- 
Bashforth predictor corrector scheme which is initialized with a fourth order Runge-Kutta 
algorithm. After propagating an arbitrarily small radiation field from the entrance to the 
exit of the resonator, a portion of the wave (depending on the reflectivity of the cavity 
mirrors) is transmitted out. of the system and the remaining wave is used to initialize the 
subsequent pass. In this manner we have been able to simulate the small signal exponential 
gain to the nonlinear saturated evolution of a tapered wiggler oscillator. Representative 
results are presented in Sec. IV. 


III. Theoretical Analysis 

The simulations are aided by monitoring a diagnostic quantity which may be obtained 
from Poynting’s theorem. (Ref. 6) 

_ Pr^)-Pr(O) .. ,_,._< 7 ( 0 )>-(*>(--» 

Vr{ : ) = - n -= 1 1e( = ) = —-:- 1 (<) 

Pb 7o - 1 


which expresses the fact, that the energy loss by the particles must appear in the radiation 
field. These quantities are monitored at. every axial grid of the simulation and provides an 
excellent check on the consistency of the field and particle evolution. 








Defining the saturated state as the condition that the radiation power gained per pass 
equals the radiation power loss at the mirrors, one obtains 


| a.(L) | 2 - | a*(0) | 2 ] = (1 - *)^™ 2 | a t (L) | 2 
Air Air 


PR»at(I J ) 


VR(L)P h 

(i -R)' 


( 8 ) 


where PR ta t is the radiation power at saturation, Pi, is the input electron beam power, rj R 

■i 

is the single pass radiation efficiency as defined in Poynting’s theorem and R, is the power 
reflection coefficient for the cavity mirrors. This expression aptly indicates the possibility 
of achieving high radiation powers in high Q resonators with modest beam powers. The 
remaining theoretical sections are devoted to estimating, t]r{L), the single pass radiation 
efficiency in terms of the prescribed beam and wiggler parameters. 

The evolution of the radiation power can be obtained from the vector potential equa¬ 
tion (4), 


d | d, 

dz 


+ »: I a. 1=^ 


*b\2fdzoFdi 


k. 



(9) 


where <f> t is the radiation phase, a $ =\ a, \ e and 0, — t/’ + 4>,- Multiplying the real part 
of Eq. (9) by | a, | (c/Air)klnivl one obtains 


dP R (z) 

dz 


Pb 


lo - 1 


kfb w | 



( 10 ) 


III. a) Resonant Particle Approximation 


The simplest estimate of the single pass efficiency as considered by Refs. 1.2 can be 
obtained by making the Compton and resonant phase approximations. The resonant phase 
approximation of Eq. (10) is given by, 


dP R (z) 

ds 


Pb 


7o ~ 1 


■ kgd, 


sin il'»R 
fit o7R 
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The sine of the resonant phase is defined by evaluating the fixed point of the pendulum 
equation (Eq. (6)) at the resonant energy given by -y R - k a fi 2 //3 z (l + 0-)k u< , with /r = 
1 4- 


n ill , 1 > d K> (14" ftz)ku<ks • I - I • , / 10 \ 

0 “ *• + *- - ~ Phi Qw 1 a ' 1 smV ' 4R - (12) 

Neglecting the rate of change of the radiation phase as compared to the wiggler period 
and inserting the results from Eq. (12) into Eq. (11) yields, 


— 77 ^ = PbVo(z), (13) 

where r]' 0 = ^jf]. Equation (13) can be integrated analytically and 

yields the following efficiency for the resonant particle approximation, 


VRet = I Vo(z)dz = 

J o 


P r (L) P r ( 0) 
Pb 


Wlo 7r(°) - 1r(L) 
(1 4- iho ) 7o — 1 


This efficiency is the relativistic generalization of the bucket efficiency of Refs. 1,2 and 
yields efficiencies which are of the order of 70 percent larger than the efficiencies obtained 
from the simulation of continuous electron beams. Although, when the simulations are 
performed for prebunched beams characterized by macro-particles injected at. the resonant 
phase, Eq. (14) yields very good agreement, with that simulation efficiency. 

III. b) Deeply Trapped Particle Compton Limit 

In order to better approximate the ensemble average over initial phases, the assump¬ 
tion of a macro-particle at the resonant phase will be replaced with an explicit average 
over the orbits of deeply trapped particles. As a first example we shall consider again the 
Compton approximation. Expanding the pendulum equation, Eq. ( 6 ), about the resonant 
phase, = 4'tR + H' at the resonant energy y r = /i 2 fc,//J ; (l 4- )k w yields, 


~ ~ - - ^1" -^ l I { sin + ros </’, R *</’}, (l- r >) 











where we shall explicitly retain the drift and acceleration of the resonant phase, rp' aR / 
0, / 0. Using the definition of the sine of the resonant phase, Eq. (12), the pendulum 

equation becomes, 


W" = -4\r{z) - *;„ n (-Ws (16) 

where k] yn (z) = (1 + /3 z )k w k a d w | a a | cos The oscillation of the deeply trapped 

particle can be explicitly solved by multiple scale techniques which exploit the desperate 
time scales of the synchrotron bouncing of the particle and the adiabatic variation of 
wiggler parameters. The solution is given by, 

~ ~t^r\ + ^o,s)cos / k ayn (z')dz'- p(z) sin / k ayn (z')dz' (17) 
'Synt-J JO Jo 

where flV’o,z) = y/k ayn (0 )(^„ - V’.fl(O) + t aR (0)/k; yn (0))/^/k ay „(z), p(z) = 

V’»/i(0)/\/^»yn(0)^»yn(“) and we have chosen boundary conditions corresponding to 
V’»(0) = V’o and V’1(0) = 0 . Making use of this solution for the ensemble average yields, 


€ -«fc/ n fc.y-t* )<<* e'L )d: Jc(p(z)), (18) 


where J m (x) is the Bessel function of the first kind of order m. By neglecting the oscilla¬ 
tions at the sideband harmonics one obtains, 


(sin^ 4 ) = sin(4’ aR (z) - 4'" R {:)/k; yn {z))J 0 (p(z))X< (19) 

with ~«/o(e(0 o,~)), 

J rfV7 0 (V) = A7o(A) + y[/r 0 (A)7 1 (A)-lir 1 (A)yo(A)], 


and 







where Hk{%) is the Struve function of order k. The results of the ensemble average, when 
included in the power evolution equation, yields, 


—= ViPb - v'i iP r( ~ )i 


where 7]) = HJ Q (p(z))cos %(- J 


I / _ 7 b (V,, Ml ■ V’,ft(~) 

and rjn = - ~^J 0 (p(r)) sin - — cos 

7o - i «; y n(') 


, , ( to ~ 1 \ 

4'sr ( = )—— -TTr 

nk V /- FF / 


Although Eq. (20) can be solved analytically for the case of constant i)' t and the 
transcendental solution for Pr(z) makes it impossible to evaluate the single pass efficiency, 
i]r(L) - ( F n(0) Pr(L))/P b- For this reason we choose to retain the slow variation of 

r)'j and i]'ji and solve Eq. (20) numerically and compare the results to the simulations. 
This efficiency yields results that, are within 10 percent of the simulation efficiency for the 
Compton regime operation. For parameters of physical interest. ij' f •*> i)' n and the evolution 
of Eq. (20) is dominated by the first term. The deeply trapped particle approximation to 
the efficiency is dominated by the bucket efficiency times the product of three form factors, 
3 , J o(Pi = ))and cos V’" R (r )/k 2 fyn (s), which are due to the average motion in the bucket , drift 
and acceleration of the bucket respectively. 

III. c) Deeply Trapped Particle Raman Limit. 

The space charge effects can be easily included with this technique of deeply trapped 
orbit, averaging. The pendulum equation in the presence of space charge waves and ex¬ 
panded about the resonant phase and energy is given by. 


- -+"r + k - -j* ~ {k ' s + 1 sin( *** 1 > 

_ sin «/-,!, - + ^, l .)l^i|c <> s(fc’, f? I p cV 

ft* ^ b r 


, (1 "b (it )k u >k t , |, | ! c , 

+- -zr— -fl.rKI COS |/', R I'll'. 

PzTr 


C*iV 


r* 

I; 
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ft: 

S’ 


ft: 

•ft 
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Sffi 
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where ip is the phase of the electrostatic potential, 4> x — and we have neglected 

the rate of change of the space charge potential as compared to the ponderomotive wave¬ 
length. By making use of the resonant phase definition in the Compton regime, k' u , - 
ks/(2fillR)da~ w /dz = (1 + /? ; )fc u ,fc,d u ,|a 4 | sin the pendulum equation can be 

cast in the form, 


H" = «(--) - P *</-, 


( 22 ) 


where a(r) = — 4'"r 


^ ^ : (A', f A- U ,)|^i| sin(V\,/? \ p 4>s 




Kyn( = ) 


2M1 zj&l 
Phi 


(k, + k u ,)\<f> x | cos(i /', R f <p - <f > s ) + 


1 + tfz )k, r k x 


n„.1 0,1 cos 


The magnitude of the electrostatic potential. |<£]|, and the phase shift between the 
electrostatic and radiation phases, y? - </>, can be estimated from the evolution equation 
for the space charge potential. 


I <t>\ 1= rr , , ( cos (V*. + V 7 - <A.<) » *>»> (</’.« + y> </>,))• (23) 

Since the left-hand side of Eq. (23) is real, one has that (sin( t/’, -f y? -</>,,)) -- 0. Using 
the deeply trapped orbits of the Compton approximation to evaluate the ensemble average 
yields, 










k ‘ ,niz) - ’k j 7 »‘\k.+k w ) C ^ J ° Ms))cos 


, (1 +fiz)k t k w . , 

+ - m ~2 -®u» °« COS V’,R, 

PllR 


where the space charge potential modifies the synchrotron wave vector k, yn by increas¬ 
ing the synchrotron period due to electrostatic repulsion, and also increases the effective 
bucket acceleration a(z). Again making use of multiple scale techniques the solution to 
the pendulum equation is given by, 


** = j^-+i(4' 9 ,z)co S f k ayn (z')dz'-p(z) S m / k ayn (z')dz', (27) 
“sj/ni 2 ) JO JO 


where £(V’o,=) = \J<k ayn (0)(if' 0 - 4'»r(0 ) - a(0)/fc^ n (0))/y fc«j,»(s) and p(z) = 

4'' a R(0)/y/i *yn(0)fc»y„(2). Performing the ensemble average with these orbits, the power 
evolution equation becomes, 


dP ^ > = ’pV'iK + s P v'nVP>>PR(-h 


where Jp »?/ = 3J 0 (/3(c)) cos = ) 

3 unl*') 


, # HR f ... v, . «(-) , / 70 - 1 V 

and ap i) n ~ •-— ~^J 0 (/?(;)) sin ---- - cos 0 jR —— ( —— ) 

l- 1 k;(z) n h V &F ) 


and 3 = J ~^J 0 (£W'o,:)) • 


The numerical solution to Eq. (28) is compared to the simulation results in the space charge 
dominated regime. The comparison presented in Sec. IV is characterized by agreement 
between theory and simulation on the order of 10 percent. 


IV. Results and Summary 


The simulation model is capable of tracking the evolution of the radiation field from 
the small signal exponential gain through the nonlinear saturated state of the oscillator. 


s. 




y v vvViV 








During the small signal evolution of the fields for which the growth rate is constant, the 
growth rate from the simulation is compared to the growth rate obtained from the solution 
of the linear dispersion relation. This comparison is presented in Fig. 1 for the case of, 
beam current I = 300.4, beam energy 7 = 2.5, wiggler field strength B w = 1000G', wiggler 
period ( w = 4.0cm, filling factor F = 0.2 and no taper of the wiggler amplitude or period. 
The comparison of the simulation and theory yields excellent agreement throughout the 
unstable spectrum. 

An example of the small signal to nonlinear evolution of the fields is presented in Fig. 
2. In this example the field evolves over 200 radiation bounce times for the case of beam 
current I = 300.4, beam energy 7 — 2.5, wiggler field strength B w = 1000G', initial wiggler 
period i w = 4.0cm, filling factor F = 0.2, linear taper of the wiggler period of 50 percent 
and mirror reflectivity R = 0.98. The saturated power and efficiency are 1 .9GIF and 17 
percent respectively; and exhibits the scaling indicated in Eq. (28). 

The single pass efficiency exhibits a weak dependence on the radiation power at the 
entrance of the resonator. In Fig. 3 we present a comparison of the single pass efficiency 
from the simulations and theory as a function of the input power. As noted in Sec. 
III. b), the theoretical evaluation of the efficiency neglects oscillations at the harmonics 
of the synchrotron period. The simulation model contains the complete dynamics of the 
electrons as dictated by the pendulum equation and the simulation efficiency tliusly exhibit s 
oscillations at synchrotron harmonics with an amplitude of the order of the well depth. 
The bars on the data points from the simulation indicate the amplitude of this oscillation 
in the efficiency. The resulting comparison indicates agreement on the order of 10 percent 
between the simulation and theoretical expressions. 

Acknowledgment. This work is supported by the U.S. ARMY STRATEGIC DEFENSE 
COMMAND. 
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